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This work describes the process of developing the nominal static aerodynamic coeffi- 
cients and associated uncertainties for the Orion Crew Module for Mach 8 and below. The 
database was developed from wind tunnel test data and computational simulations of the 
smooth Crew Module geometry, with no asymmetries or protuberances. The database 
covers the full range of Reynolds numbers seen in both entry and ascent abort scenarios. 
The basic uncertainties were developed as functions of Mach number and total angle of 
attack from variations in the primary data as well as computations at lower Reynolds num- 
bers, on the baseline geometry, and using different flow solvers. The resulting aerodynamic 
database represents the Crew Exploration Vehicle Aerosciences Project’s best estimate of 
the nominal aerodynamics for the current Crew Module vehicle. 
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c x 

A force or moment coefficient 

C A 

Axial force coefficient 

Cd 

Drag force coefficient 

C l 

Lift force coefficient 

Cm 

Pitching moment coefficient 

Cm 

1 1 1 apex 

C m , resolved at vehicle theoretical apex 

Cm cg 

C m , resolved at cg location 

C m 

ni cqx 

C m , resolved at mrc along centerline 

Cn 

Normal force coefficient 

C n 

Yawing moment coefficient 

eg 

Center of gravity 

d 2 

Bias correction factor for range-based 
standard deviation 

FMV 

Blended velocity function 

k 

Coverage factor, y/3 

kfps 

Thousand feet per second 

L/D 

Lift-to-Drag ratio 

M 

Database cardinal Mach number 

M 

Margin Index 


Mqo Freestream Mach number 
mrc Moment reference center 

Rcb Reynolds Number based on vehicle 

diameter 

Ri Range of data for condition i 
sf Smoothing factor 

uC x Uncertainty in a force or moment coefficient 

V Fc x Database uncertainty factor for a force or 

moment coefficient 
Uoq Freestream velocity 

a Angle of attack, deg. 

ar Total Angle of attack, deg. 

/3 Sideslip Angle, deg. 

cj) Roll Angle, deg. 

S Increment or range to cover 

a Standard deviation 

Subscripts 

5sp Five species air 

alt Altitude 
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bal 

Balance accuracy 

Itp 

Laminar-turbulent variations 

c2c 

Overflow to Usm3d comparisons 

ovf 

Overflow 

cfd 

CFD variations 

P9 

Perfect gas air 

chem 

Chemistry model differences 

Re 

Reynolds number variation 

fit 

Flight Reynolds number 

RR 

Repeatability 

hs 

Heatshield asymmetry 

turb 

Turbulent 

idat 

IDAT backshell angle change 

wt 

Wind tunnel data variations 

interp 

Interpolation 

wt2cfd 

Wind tunnel to CFD variations 

lam 

Laminar 

wtRe 

WT Reynolds number 


I. Introduction 


SM 

Service Module 

CM 

Crew Module 


SA 

Spacecraft Adapter 



LAT 

Launch Abort Tower 


Figure 1: Orion Crew Exploration Vehicle components. 


The Apollo-derived Orion Crew Exploration 
Vehicle (CEV) was designed by NASA and its 
industry partners within the now-cancelled Con- 
stellation Program as part of the Agency’s Ex- 
ploration Mission, and was intended to be the 
foundation for manned exploration of the Moon, 

Mars, and beyond. 1,2 The Orion CEV design is 
now the reference vehicle for the development of 
the Multi-Purpose Crew Vehicle (MPCV), the ex- 
ploration vehicle that will carry crew to space, 
provide emergency launch abort capability, sus- 
tain the crew during the space travel, and provide 
safe re-entry from deep space return velocities. 3,4 
The CEV (and now the MPCV) consists of the 
Crew Module (CM), Service Module, Spacecraft 
Adapter, and Launch Abort Tower, as shown in 

Figure [l] The Orion CM is similar in shape to, but larger than the Apollo capsule. A test flight planned for 
late 2013, designated Orion Flight Test 1 (OFT-1), will focus on the entry phase of flight for the CM. 

The Orion aerodynamic database 5,6 has been developed 
by the CEV Aerosciences Project (CAP), and is regularly 
updated with improvements to the aerodynamic modeling of 
various systems. The primary function of the database is to 
provide aerodynamic data to the trajectory simulations that 
are used to develop the guidance, navigation, and control 
systems for the vehicle and provide targeting and landing 
ellipse prediction during flight operations. The aerodynamic 
database development process is shown notionally in Fig- 
ure [2] Note that this paper will use the term aerodatabase 
throughout to refer to the Orion aerodynamic database. The 
CAP team provides this data through an API (Application 
Programming Interface) that is integrated into the trajec- 
tory simulation tools. The API uses tabulated nominal and 
uncertainty aerodynamic data to compute and return the 
aerodynamic forces and moments acting on the vehicle at 
the desired vehicle state. The tabular data is developed from 
various computational and experimental sources. Uncertain- 
ties due to turbulence modeling, grid resolution, wind tunnel 
repeatability, and other physical modeling are combined to 

provide tabulated database uncertainties. The CAP database covers the aerodynamics for all phases of the 
vehicle flight beginning with the separation from the launch system (including nominal and abort situations) 
until the CM re-enters the atmosphere and lands. 

This paper describes the general process of developing the nominal static aerodynamic coefficients and 
uncertainties for the subsonic through low hypersonic flight regimes for the Orion Crew Module (M^ < 8), 




Figure 2: CAP aerodatabase development pro- 
cess. 
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The process described reflects the current best estimate of the database, and discusses areas where further 
work is planned to support the first entry test flight of the Orion vehicle. This paper covers the full angle- 
of-attack range for the vehicle, 0° < a < 180°, for < 8. A companion paper 7 covers the database 
development for > 8. 

The available data, both experimental (WT) and computational (CFD), are discussed, with particular 
attention paid to comparisons between the various data sets. The basic process for combining the available 
data into a single data set varying in Mach and angle of attack is outlined, and some examples of the 
blending will be shown. The formulation of the uncertainties and development of each term is discussed, 
with particular attention paid to the Reynolds number variations and CFD to WT comparisons. The 
final database nominal coefficients with uncertainties are compared to the available data and vehicle trim 
characteristics are presented. 



Figure 3: Dimensions of the smooth, axi-symmetric 
baseline Orion Crew Module, in inches. 


II. Database Inputs 

II. A. Orion Crew Module Geometry 

The nominal analytical Orion CM geometry is based on 
the Apollo configuration, and is shown in Figure [ 3 ] The 
spherical heatshield and conical backshell have been 
scaled to a maximum diameter of 198 in compared to 
Apollo’s 154 in. The CEV apex is truncated to accom- 
modate docking hardware. 

The flight geometry is still being developed, and de- 
parts from the nominal, axisymmetric geometry in sev- 
eral key areas. The aerodatabase addresses these geom- 
etry differences by incorporating additional analysis to 
adjust the nominal coefficients and using uncertainties 
to cover expected variations. The geometry variations 
fit into three main categories. 

First, the nominal 32.5° backshell angle was 
widened by 2.5° (to 30°), moving the theoretical apex 
further from the vehicle base, as shown in Figure [4j 
This modification provides more packaging volume for 

the parachute system and is referred to as the IDAT geometr^] Initial experimental and computational 
data were generated on the axisymmetric nominal geometry. Subsequent computational studies on the IDAT 
geometry have been performed, and the effect of incorporating the backshell angle change into the database 
is covered below. 

The second major departure from the nominal geometry was reshap- 
ing the heatshield to minimize the thickness of the thermal protection sys- 
tem (TPS), resulting in an asymmetric heatshield shape. Conceptually, the 
shape is designed to be thicker in the higher heating regions such that the 
expected ablation drives the shape closer to the nominal spherical shape, 
and implies that the effect of the asymmetry will be greatest at the high- 
est entry speeds. While this shape is still evolving, there have been some 
studies to address the aerodynamic effect. The studies have been primarily 
at hypersonic conditions, and are briefly discussed below in the uncertainty 
section. 

The last group of geometry variations include features such as footwells, 
windows, steps in backshell tile thicknesses, and other protuberances. The 
aerodynamic effects of these relatively small features have not been quanti- 
fied, and are assumed to be accounted for within the uncertainty model, as 
discussed further in the uncertainty section. 



Figure 4: Dimensions of the 

axi-symmetric baseline Orion 
Crew Module. 


a The effort to redesign the packaging and deployment schedule for the parachute system was called the Integrated Design 
Assessment Team, and the resulting changes to the CM vehicle are generically called the IDAT geometry. 
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II. B. Crew Module Coordinate System Conventions 

The CM coordinate system, angle of attack, and aerodynamic coefficient orientation conventions are shown 
in Figure 0 taken from the CEV Aerodynamic Databook. 5 Note in particular that a heatshield-forward 
attitude has an angle of attack of 180°. 


Z-axis 



Figure 5: Axis, Force, and Moment Definitions for Crew Module . 5 Astronaut, shown for orientation purposes, is not to 
scale. 


The aerodynamic moments are typically resolved about a nominal eg location defined by the Orion 
Project. Through the course of the CM development, the nominal eg location has shifted, primarily providing 
progressively smaller offsets of the eg location in the z-axis direction. The CAP team typically considers 
moments resolved about three different mrc locations: apex, eg , and symmetric eg. The database is provided 
with the moment about the theoretical apex for the nominal, 32.5° backshell, geometry (labeled MRC in 
Figure 0). The flight eg location is used to determine flight characteristics such as trim angle of attack, and 
is provided as an output of the API. The symmetric eg (a location along the x-axis corresponding to the eg 
with no y- or z-axis offsets) is used to develop the nominal pithching-moment coefficient so that symmetry 
conditions can be enforced. Moment uncertainties are developed for either the eg or symmetric eg location. 


II. C. CM Database Formulation 

The CM portion of the database provides the aerodynamic forces and moments as a function of a velocity 
parameter, FMF, and the orientation of the vehicle with respect to the flow. The database formulation 
takes advantage of the fact that the vehicle is primarily axisymmetric and treats the nominal coefficients as 
functions of FMV and olt only. Tables are provided for Ca, Cjv, and C m , and transformations made to 
compute the full set of coefficients. The velocity parameter, FMV , is defined in Equation [l] as 


Moo 

FMV = { Moo(9.8 - Uoo) + UooiUoo - 8.8) 

Uoo 


if U oo < 8.8 kfps 

if 8.8 kfps < Uoo < 9.8 kfps 

if Uoo > 9.8 kfps , 


(i) 


where Uoo is given in units of kfps. The FMV function provides a single velocity parameter for the entire 
database by taking advantage of the near numerical equivalency of Mach number and velocity (in kfps) around 
Mach 1(0 The blended parameter provides a smooth transition between Mach number for subsonic and 
supersonic speeds and velocity for hypersonic speeds. For the development of the database for FMV < 8, 
Mach number and FMV are used interchangeably. 

b Note that the formulation of FMV complicates the direct application of the Orion aerodatabase to entry simulations in 
other atmospheres, such as at Mars. 
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Because the vehicle is predominately axi- symmetric, the database formulation treats the nominal coef- 
ficients as functions of FMV and olt only. Tables are provided for Ca, CV, and C m , and transformations 
made to compute the full set of nominal coefficients. 

The database API computes dispersed aerodynamic coefficients based on the data in the uncertainty 
tables. In order to facilitate development of dispersed trajectory simulations (typically Monte Carlo based), 
the user provides an uncertainty factor for each uncertainty coefficient, and the nominal and uncertainty val- 
ues are combined to form the dispersed coefficient. For a particular simulation within a dispersed trajectory 
simulation, the user provides an uncertainty factor, typically in the range — 1.0 < UFc x < 1.0, which will 
be applied to the aerodynamics for that simulation. An uncertainty factor of zero will return the nominal 
coefficient for C x , and UFc x = 1.0 will return the nominal plus the total uncertainty for C x . For the Orion 
database, all uncertainties are specified as uniform uncertainties except for the rolling moment uncertainty 
which is treated as a normal distribution. This means that the uncertainty factors chosen for the dispersed 
trajectory set will be chosen based on a uniform distribution from —1.0 to +1.0. 

Table [l] specifies the database tables required for the CM static aerodynamics, and lists the independent 
parameters for each table. Details of how the data from the tables are built into the final aerodynamic 
coefficients can be found in the companion paper, 7 and more fully in the Orion Formulation Document 6 

Table 1: Database tables and arguments for the CM static aerodynamics. 6 



Nominal Tables 


Uncertainty Tables 


Coeffi- 

cient 

Table 

Name 

Table 

Arguments 

Coeffi- 

cient 

Table 

Name 

Table 

Arguments 

C A 

CACM 

FMV qlt 

uCd 

UCDCM 

FMV 

OLx 

Cn 

CNCM 

FMV olt 

uC l 

UCLCM 

FMV 

OLx 




u(L/D) 

UL0DCM 

FMV 

OLt 

C m 

• ' 0 apex 

CMCM 

FMV c^t 

u C mcg 

UCMCGCM 

FMV 

a /3 




uC„ cg 

UCLNCGCM 

FMV 

a (3 




uCi 

L cgx 

UCLLCM 

FMV 



II. D. Available Experimental Data 

The primary experimental data available for M ^ < 8 come from three tests conducted early in the project: 
05-CA 8 in the Ames Unitary Plan Wind Tunnel (AUPWT), 03-CA 9 in the Langley Unitary Plan Wind 
Tunnel (LaUPWT), and 09-CA 8 in the Langley Aerothermal Laboratory’s Mach 6 facility. There have been 
several additional tests, for the LAV configuration and dynamic CM testing, that have generated static 
aerodynamic data on the CM, and have been used in developing uncertainties. A summary of the available 
test data on the CM is given in Table [2] 

The 05-CA test used both a small (3.0%-scale) and large (7.7%-scale) model tested in both the the lift 
and 9x7 legs of the AUPWT. The test covered an angle of attack range from 142° to 172°, for 0.3 < < 2.5 

and l.OxlO 6 < Rejj < 5.3xl0 6 . Additional data was obtained for M ^ = 0.5 at Rep = 7.0xl0 6 , 7.6xl0 6 and 
for Moo = 1.6, 2.5 at Rep = 0.5xl0 6 . The data from 05-CA with the large model at the highest available 
Reynolds number were used in the database development. The 03-CA test utilized the 3.0%-scale model 
tested in the LaUPWT over an angle-of-attack range of 140° to 170° for Mach numbers 1.6, 1.8, 2.0, 2.5, 
3.0, and 4.5, and Reynolds numbers of 0.5xl0 6 , l.OxlO 6 , and 1.5xl0 6 based on body diameter. The data for 
Moo >2.5 and Rcd = 1.5xl0 6 from 03-CA are utilized in the database. The only hypersonic aerodynamic 
test, 09-CA, was run in Langley’s 20-Inch Mach 6 Air facility, and covered the full range of angle-of-attack, 
-5° to 185°, and Re d = l.OxlO 6 . This test was run with an apex cover on to provide a geometry similar to 
Apollo, and without the apex cover for the CEV configuration. Murphy et al 10 cover the details of these 
tests and provides discussion of the data including Mach and Reynolds number effects that guide how the 
data is used in the database. 
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Table 2: Available experimental data for CEV CM, Moo < 8.0. 


Test 

Ref 

Facility 

Mach 

a range 

Ren, 10 -6 

notes 

03-CA * 

9,10 

LaRC UPWT 

1.6 - 4.0 

140°-170° 

1.0 - 3.0 

3.0% scale 

05-CA * 

8,10 

Ames 11 ft 
Ames 11 ft 
Ames 9x7 

0.2 - 1.4 
0.5 

1.6 - 2.5 

142°-172° 

142°-172° 

142°-172° 

max 5.3 
max 7-7.5 
max 5.3 

7.7% scale, 3% scale 
7.7% scale, 3% scale 
7.7% scale, 3% scale 

09-CA * 

10,11 

LaRC 20-Inch 
Mach 6 

6.0 

-5°-185° 

1.0 

2.02% scale, apex 
cover off / on 

24- AA 

12 

AEDC 16T 

0.3 - 1.2 

142°-172° 

1.0 

Reo too low, asymmet- 
ric CM geometry 

25- AA 

12 

Ames 

1.6 - 2.5 

142°-172° 

0.75 

Ren too low, asymmet- 
ric CM geometry 

27-AD 

13,14 

LaRC TDT 

0.2 - 1.1 

142°-172° 

1. 0-5.0 

Dynamics test, sting 
interference issues 

61-AA 

15 

LaRC 14x22 

0.13 

142°-172° 

1.0 

Reo too low, asymmet- 
ric CM geometry 


Data used for CM static aerodynamic database development. 


II. E. Available Computational Data 

The available computational data for the CM comes from several CFD solvers and covers both the baseline 
geometry with the 32.5° backshell angle and the ID AT geometry with the 30° backshell angle. The data sets 
with their respective Mach number and angle- of- attack ranges are summarized in Table [ 3 ] 

Table 3: Available computational data for CEV CM, M^ < 8.0 


CFD Solver Backshell Angle Mach number range a range 


Overflow 
Overflow 
Usm3d 
Laura, Dplr 


32.5° 

30° (IDAT) 
30° (IDAT) 
32.5°, 30° 


0.2, 0.3, 0.5, 0.7, 0.9, 1.1 0°- 180° 

0.3 - 8.0 0°- 180° 

0.5, 0.7, 0.9, 1.2, 2.5, 8.0 0°- 180° 

2.0 - 8.0 150°- 172° 


There are 2 primary sets of CFD data for the CM, both generated with the Overflow 16-19 solver using 
the SST turbulence model. Both primary data sets cover the full angle of attack range, 0° < a < 180°, with 
the point distribution highest in the trim region. The first set is for the nominal axisymmetric geometry, and 
are available for M 0 0 = 0.2, 0.3, 0.5, 0.7, 0.9, and 1.1. There is a corresponding set of computations on the 
IDAT geometry with the wider backshell angle. These computations cover the full Mach range to = 8.0. 
The IDAT geometry cases were used for the nominal coefficient development. 

The nominal backshell angle data set includes computations at both wind tunnel and flight Reynolds 
numbers, as well as some with a sting configuration matching the 27- AD test. There was considerable 
duplication in this set, with cases run on similar grids with the same solver options by multiple researchers. 
These duplicate cases were treated as repeat solutions to aid in the development of uncertainty. The cases 
with multiple Reynolds numbers were also used in the uncertainty development to cover the range of Reynolds 
numbers seen in flight, particularly subsonically. 

An additional set of computations at selected Mach numbers and angle-of-attack ranges for the IDAT 
geometry were made using the unstructured package TetrUSS, 20-22 which is comprised of the Vgrid 
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meshing package and Usm3d flow solver. The flow conditions overlap with the Overflow IDAT set and 
the computations have been utilized in the uncertainty development. 

As Mach number increases into the hypersonic regime, real gas effects become a significant factor in the 
aerodynamics; these are not modeled either in the wind tunnel or with the primary Overflow data sets. 
For the higher Mach numbers (M^ > 2), a limited set of computations from the Laura 23,24 and Dplr 25 27 
solvers has been used for nominal coefficient development in the trim region where appropriate, and for 
uncertainty quantification for real gas, altitude, and laminar-turbulent variations for 2 < M ^ < 8. This 
data set was developed by the CAP Aerothermal team as part of their efforts to characterize heating on 
the CEV, and is referred to in this paper as the aerothermal data set. These computations cover velocities 
above Mach 2, angles of attack between 150° < a < 172°, and a range of altitude conditions. For M ^ > 8, 
the aerothermal data is the primary data source, and so including this data at the higher Mach numbers 
facilitates blending between the two regions of the aerodatabase. Note that the data used for nominal and 
uncertainty development is the average of a laminar and a turbulent solution at the same conditions, and is 
referred to in this paper as a Itp data point. More details on the full aerothermal data set can be found in 
the companion paper 7 which covers the hypersonic aerodynamic database development. 

II.F. Data Coverage, Angle of Attack 

The aerodatabase by design covers the entire range of orientation for the CM, such that a varies from 0° to 
360°, and /? from -90° to 90°. When the coefficients are axisymmetric, cut varies between 0° and 180°. 

Figure [6(a)] shows the various available data sets, computational and experimental, for the full range of 
Mach number and angle-of- attack. The trim region for < 1.6 is shown in Figure [6 (b)| The gray symbols 
and lines show the points in the database tables for the portion of the CM database covered in this paper. 
The wind tunnel data are shown in thick lines in shades of blue. The Mach 6 data from the 09-CA test 
can be seen to extend over the full angle- of- attack range, while the 05-CA and 03-CA data are only in the 
trim region, for < 4. The CFD for the baseline geometry, shown by the blue squares, are duplicated in 
most, but not all, cases by the IDAT CFD data. The orange diamonds are the primary Overflow IDAT 
CFD data used for the nominal development. The purple triangles are the duplicate points computed with 
Usm3d, and are only available for selected Mach numbers and limited angles of attack. 


■ ■■ ■■■ ■MHMHM g mmm 


>000 0 000 




00000004 ♦ ♦♦♦ 
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► 444 4 4 4 4 43404010004 4 
44444444 M 44444444440 
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Angle of Attack, deg. 
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A CFD, USM3D, IDAT 
O Aerothermal CFD 
Database 



(a) Moo <9, 0° < a < 180° (b) M^ < 1.6, 130° < a < 180° 

Figure 6: Mach number and angle-of-attack range coverage for all available data. 


Typical nominal entry and abort trajectories fly in a tight angle-of-attack range, trimming in a 10° range 
subsonically and a 7° range hypersonically with full uncertainties applied. This range is located somewhere 
between 145° and 175° depending on the nominal eg location, with the nominal trim angle of attack closer to 
180° as the z-cg location is closer to the vehicle centerline. The clustering of CFD points and concentration 
of the experimental data reflect the need for the best definition of the aerodynamics for this angle-of-attack 
range. At the lower subsonic speeds, the dynamics become very important, with the vehicle experiencing 
large oscillations about the trim point. The total angle of attack, cut, can be as low as 120° below Mach 0.5. 
This implies that there is a broad range of angles-of-attack and sideslip angles that must be well-covered in 
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the lower speed ranges, and therefore a denser CFD angle-of-attack distribution for subsonic Mach numbers. 
Below a = 45° there is a very sparse distribution in CFD data. For normal operations, the vehicle should 
not fly in this region. 


II. G. Reynolds Number Coverage, Trajectory Envelope 

The available data is limited in the variation in Reynolds number, particularly for the experimental data. 
Historical Apollo data 28-36 and the current CFD computations suggest that variation in the aerodynamics 
(particularly drag) continues as Reynolds number increases beyond the Rep = 5.3 xlO 6 available with the 
CEV WT data for subsonic Mach numbers. Figure [7] shows the drag coefficient for the current data set 
at Mqo = 0.7 compared to relevant Apollo WT and flight data, for a > 140°. The primary 05-CA data, 



— Database Nominal 
0 05-CA, Re D = 1.0M, 1.5M 

0 05-CA, Re D = 5.3M 

A CFD, 32.5°, Re D =21M 

A CFD, 32.5°, Re D =4-5M 

^ CFD, IDAT, R^=21 million 

v CFD-USM3D, IDAT, R^=21M 

ApolloWT, M=0.70, Re D =1M 

+ ApolloWT, M=0.70, Re D =5M 

ApolloWT, M=0.72, Re D =8M 

+ ApolloWT, M=0.72, Re D =13M 

A AS-202 

* Apollo 6 

# Apollo 7 


Figure 7: Available experimental and computational data for — 0.7, compared to historical Apollo WT and flight 
data. 

at the highest WT Reynolds number, are shown as the cyan circles, and lie in the middle of the available 
data. The lower Reynolds number data from 05-CA, yellow circles, has a higher drag. The Apollo WT 
data (+ symbols), collected from various tests, is consistent with the current CEV data, with the higher 
Reynolds number data showing lower drag. Note that the Ren = 13xl0 6 data is mingled with the lower 
Re^ = 8xl0 6 data, and that there is a large scatter in both sets. The Apollo flight data, from 3 of the 
early test flights, shows higher drag than most of the data. This pattern is consistent with other Mach 
numbers. The CFD data has generally lower drag than the experiment, and the low Reynolds number CFD 
shows higher drag than the high Reynolds number CFD. This pattern shown here for = 0.7 is consistent 
with other subsonic Mach numbers, and begins to show the difficulties in modeling the aerodynamics in the 
angle-of-attack range in which the vehicle trims. Differences between the current WT data and the CFD are 
discussed in the next section. 

The Reynolds number range for the available computational and experimental data are shown in Figure [8j 
with representative entry and abort trajectory datsj^] shown by the dark and light gray lines, respectively. 
Because the vehicle must be capable of performing an abort at any point along the ascent trajectory, the 
range of altitude conditions that the CM might experience is large. This is shown in Figure [8] by the large 
spread in Reynolds number at any given Mach number. The nominal entry trajectories have a much smaller 
spread in Reynolds number. Note that trajectories shown here are intended to show general trends, and are 
not necessarily equivalent to the final flight profiles. 

c Monte Carlo sets of trajectories provided by the Orion Ascent Abort and Entry Mode Teams. 
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Figure 8: Reynolds Number ( Reo ) vs. Mach number for selected entry and abort trajectory sets. 


The available WT data, with a maximum Re d = 5.3 xlO 6 (except at = 0.5), is below the nominal 
entry trajectory Reynolds number for all Mach numbers, and on the bottom edge of the abort trajectories 
subsonically. This lack of coverage of the trajectory space in Reynolds number is particularly problematic 
for the subsonic Mach numbers, were the CFD data predicts significant differences between the WT and 
flight Reynolds number computations. This suggests that the WT data at the highest Reynolds numbers 
have not reached Reynolds number independence. The problem is compounded because the CFD does a 
poor job of matching the WT data for the WT conditions. The approach for the database development has 
been to use the highest Reynolds number WT data as if it were equivalent to the flight Reynolds number, 
and account for a level of Reynolds number variation within the uncertainities. 

The lack of coverage in the WT data for flight Reynolds number is less of a problem above Mach 1. The 
WT data show less variation with Reynolds number, and much of the lower range of the trajectory space is 
covered by the data. Above Mach 1, the primary CFD is available only for the flight Reynolds number. 

Above Mach 3, the aerothermal data set covers the Reynolds number range fairly well, due to the range of 
altitude conditions considered by the aerothermal team. As will be shown later in the paper, the magnitude 
of the variation in aerodynamic quantities due to altitude variation is small, but measurable, and is accounted 
for in the uncertainty model. 

II. H. Subsonic WT and CFD Data Comparisons 

For the CM, the subsonic trim region has been the most difficult aerodynamic environment to develop. As 
discussed in the preceding sections, the wind tunnel conditions cover a much lower Reynolds number range 
than is seen in flight, and this prompts concerns over how well the data extends to actual flight. The CFD 
data show large discontinuities in drag at particular M^-a combinations that are not seen in the wind tunnel 
data, and this lowers confidence in how well the CFD is modeling the large, unsteady, separated base flow 
region which is dominated by difficult to model turbulence effects. 

Figure [9] shows the available drag coefficient data for several Mach numbers. For M ^ = 0.3, the CFD 
at WT Reynolds number (left triangle) compares well with the experimental data, and is higher than the 
CFD at flight Reynolds number. There is significant variation between the repeated CFD cases (on the 
baseline geometry) as well as for the IDAT geometry. For M ^ =0.5, the CFD results are even more varied, 
and the Overflow IDAT computations show a jump in drag around a = 155°. The differences between 
the computations at WT and flight Reynolds number again show a significant difference, although the flight 
Reynolds number data line up with the WT data for this Mach number. At = 0.9, the Reynolds number 
differences are small, but the jump in drag is seen in all of the computations except those with Usm3d. The 
discontinuity in drag goes away above Mach 1. At M ^ = 1.2, the computational data is smooth again, but 
the drag is higher than the WT. As Mach number increases further, the agreement between WT and CFD 
data are improved, and the CFD data does not exhibit the instabilities seen at lower Mach numbers. 

The details of the computational solutions and the ongoing efforts to improve the agreement of the 
solutions with the WT data are explored in both Stremel et al 37 and McMillin et al 38 for the Overflow 
and Usm3d computations, respectively. Murphy et al. 9 10 provides more detail on the subsonic Reynolds 
number variations in the experimental data. At this time, a wind tunnel test in NASA Langley’s National 
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Figure 9: Drag Coefficient vs. angle of attack for available experimental and computational data for a > 130° at a range 
of Mach numbers, M oo = 0.3 to 1.6. 


Transonic Facility is planned to provide flight Reynolds number data for M ^ <0.9. 

As a result of the difficulties with the CFD solutions, the database nominals for M ^ <1.5 are developed 
where possible from the wind tunnel data only. The CFD data is used heavily in the uncertainty development 
and for the nominals outside of the 140° < M ^ < 172° range. 
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III. Nominal Coefficient Development 


The nominal coefficients for the database have been developed by smoothing and merging the avail- 
able data sets at the cardinal database Mach numbers (M ) using curve fitting techniques implemented 
in MatLab ®. 39 More advanced response surface and data fusion methodologies were not used initially in 
order to focus on the data quality issues at the subsonic Mach numbers. This simple methodology was 
then expanded to include breakpoints for M < 8. It is anticipated that the more advanced methods will be 
implemented in the next database update. 

The basic smoothing and blending techniques common to the development are described first. Then the 
nominal development discussion is grouped according to the methodologies applied to each block of Mach 
numbers. The smoothing and blending techiques have been applied and adjusted to suit the available data 
in a given database space. The experimental data dominant block covers the Mach numbers below 1.4 where 
both WT and CFD data is available. Here, the WT data is used exclusively for the range of angles attack 
where it is available, and blended to the CFD data outside of this range. In the CFD + WT combined block, 
1.6 < M < 4, the WT and CFD data are combined over the WT angle-of-attack range. There are three 
Mach numbers (M = 0.6, 0.8, 1.4) where experimental data are available that were not originally cardinal 
Mach numbers in the database, and no CFD has been run at these conditions. The WT data suggested 
that the trends in Mach number were not linear between the existing cardinal Mach numbers, and database 
breakpoints were developed in this WT + interpolated CFD block. The nominals for M = 6 are developed 
with the CFD and WT data combined for a > 140°, but only use the CFD for a < 140° so as to better 
capture the effects of the wider IDAT backshell angle. The last set of nominals is for M = 8.0. These 
combine the aerothermal CFD to the primary set of Overflow CFD data in a similar manner as for the 
1.6 < M < 4 conditions. 

III. A. Smoothing and Curve Fitting Methodologies 

The basic process for developing the nominal coefficients for the database was to develop smooth data sets 
for the experimental and CFD data separately. The two sets were then combined according to how each set 
was to be used within a particular Mach number range. 

The general smoothing approach uses the f SmoothingSpline ; type of fit from the MatLab® curve 
fitting toolbox. 39 The methodology provides one tuneable parameter, the smoothing factor (sf). This 
relaxation parameter controls the degree to which the data is forced toward individual data points and can 
vary between 0 and 1. A value of 1 should force the curve to go through all data points, while small numbers 
produce a looser fitting curve. For smoothing both the CFD and experimental data, sf was set to 0.1, unless 
otherwise noted. For combining the CFD and WT data, sf was set to 0.01 unless otherwise noted. 

The data for each of the four CFD data sets were separately smoothed in angle-of-attack into single 
curves for each for each coefficient (Cd, Cl, Cm ) at each Mach number, defined at 1° increments over the 
full angle-of-attack range. For the 05-CA and 03-CA experimental data, the smoothing approach was used 
to combine runs where repeat run data were available, and to provide a smooth curve with the desired 1° 
increment in angle-of-attack. Note that the experimental data was only available over the angle of attack 
range of 142° < a < 170°. The data from the Mach 6 test (09- C A) was handled different 1}Q The highest 
Reynolds number runs were selected for each Mach number, and points at similar angles of attack were 
grouped, The angle of attack, Cd, Cl, and Cm cgx were averaged for each group, and then Cn, Ca, Cm , C mcg 
and L/D were computed from the averaged data. Sample smoothed curves for both CFD and WT data 
are shown in Figure [lOj for = 0.5, along with the averaged data for Mach 6. The smoothed data CFD 
and WT sets were then combined to form the nominal coefficient data at each Mach number, and further 
manipulated for a < 120° to return the nominal curve to the averaged CFD data, rather than the smoothed. 
These steps are described in more detail in the next sections. 


d The 09-CA data had been merged into a single set for the database much earlier in the project, and this merged set was 
used herein rather than the smoothing approach used for the lower speed tests. 
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Figure 10: Examples of smoothing process for CFD and WT data sets. 


III.B. Experimental data dominant, M = {0.3, 0.5, 0.7, 0.9, 0.95, 1.05, 1.1, 1.2} 

For the subsonic and transonic Mach numbers, wind tunnel data was used exclusively where available because 
of the large variations seen in the CFD solutions in this region and the low drag overall predicted by the 
CFD, as discussed earlier. Before the smoothed WT and CFD data were combined for the final nominal 
curve, the CFD curve was modified to allow the WT data to dominate in the trim region. First, the data 
points in the range 142° < a < 170° were removed. Adjustments to the CFD curve were then made for the 
range of 170° < a < 180°. For Cd and Ca, the CFD curve (for a > 170°) was anchored to the WT curve 
at a = 170°, so that the trend of the CFD data was used to fill in to 180°. For Cl , CV, and C m , the CFD 
curve was retained, with the coefficient value forced to zero at a = 180° so that the final curve would honor 
symmetry conditions at a = 180°. For all coefficients, points for 171° < a < 178° and 130° < a < 142° were 
then removed, so that the smoothing process could better handle the transition between CFD to WT data. 
For a < 120°, the CFD curve described in the previous section was retained. The smoothed WT curve and 
the modified CFD curve were then combined into a final nominal curve using the MatLab® ‘smoothing 
spline’ with the sf set to 0.01. Figure [Tl] graphically shows this process. The points shown for the smoothed 
WT (blue circles) and modified CFD (red diamonds) are the data that are used to produce the combined 
curve (black). 

To finalize the nominal curve, point density of the blended curve is reduced to match the database table 
breakpoints in angle-of-attack (every 5° between 0° < a < 130° and every 2° between 130° < a < 180°) 
Additionally, the data for a < 120° was replaced by the original averaged CFD data. When there was not 
a data point at a specific database angle of attack (5°, 10°, etc.), the value from the smoothing process was 
retained. 

Figure [12] shows the original WT and CFD data, the combined curve developed from the modified CFD 
curve, and the final nominal curve for Cd and Cl for the M = 0.5 condition. The effect of both the reduced 
density in the curve definition and replacement with the original CFD data is particularly noticeable at the 
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Figure 11: Blending of WT and CFD, where WT data are dominant, M = 0.5 


lowest Mach numbers near a = 60° where there is a steep gradient the Cl (and C rric , not shown) data, 
Figure [12(b) | 




Figure 12: Final nominal curve over full angle-of-attack range for M = 0.5 


III.C. CFD + WT Combined, M = {1.6, 2.0, 2.5, 3.0, 4.0} 

For the database breakpoints from M = 1.6 to M = 4.0, the experimental and CFD data are combined 
equally, as discussed in Section m Note that the WT data was developed for the baseline geometry, and 
the CFD for the IDAT geometry. There are very limited CFD comparisons that suggest there is minimal 
difference in the aerodynamics for the baseline and IDAT geometries for > 1.6 in the WT angle-of- 
attack range, and so the influence of the WT data is maintained. The process used is similar to that for 
the WT data dominant described in the previous section, except that the smoothed CFD data is retained 
over the angle-of-attack range where there is WT data. The CFD curve is still modified by deleting data 
between 130° < a < 142° to facilitate blending. The CFD curve is again anchored (for CU, Cl > ) at 170°, this 
time to the average of the WT and CFD data, and then the points for 171° < a < 178° removed from the 
CFD curve. The modified CFD curve and the smoothed WT data are then combined using the MatLab® 
smoothing process to provide the final nominal curve. As with the previous process, the CFD smoothed 
curve is replaced with the averaged data for a < 120° when the curve is coarsened for the final database 
curve. Figure [T3| graphically shows this process. 
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Figure 13: Blending of WT and CFD, where WT and CFD data are combined equally, M = 2.5 


III.D. WT + Interpolated CFD, M = 0.6, 0.8, 1.4 

In versions of the aerodat abase prior to the most recent, the CM database tables included Mach number 
breakpoints at M = {0.3, 0.5, 0.7, 0.9, 0.95, 1.05, 1.1, 1.2, 1.6}, and CFD cases were typically run for these 
Mach numbers. The 05- C A WT data includes runs at = 0.6, 0.8, 1.4 which were not previously incor- 
porated into the database. Linear interpolation is used by the database between breakpoints, but the trends 
in the WT aerodynamic coefficients are not linear with Mach number. Figure [14] shows the variation in 
drag coefficient with Mach number for a range of angles of attack. The non-linearities are greatest between 
Moo = 0.5 and =0.7 for drag, and strongest between =0.5 and M^ = 0.7 for pitching moment^ 




Figure 14: Cd vs. Mach number for several angles of attack, smoothed 05-CA wind tunnel data. 

To include these additional Mach number breakpoints in the database required creating a CFD curve by 
linearly interpolating in Mach number for each angle of attack between the surrounding available CFD curves, 
as shown in Figure \lE\ for M = 1.6. The interpolated CFD curve was then combined with the WT data in the 
same manner as was done for the M < 1.2 breakpoints, with the CFD ignored between 130° < a < 178°. 
Note that outside of the region where the wind tunnel data were available, no attempt at using the trends 
from the trim region was made; the current approach yields similar results as the linear interpolation used 
by the database when the additional breakpoints are not included. Figure [T6| shows the final database curve 
for M = 0.6 plotted with both the WT data at = 0.6 and the CFD data for = 0.5, 0.7 for a > 110°. 


e Murphy et. al. 10 includes thorough discussion of Mach effects. 
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Figure 15: Examples of combining CFD at nearby Mach numbers to get M = 0.6 CFD curve 




Figure 16: Examples of combining interpolated CFD with WT to get final M = 0.6 nominal curve. 


III.E. Hypersonic and Blending with WT, Mach 6 

Both WT and CFD are available for the full angle-of- attack range for M = 6. However, since the wider 
backshell angle for the IDAT geometry would be expected to influence the aerodynamics for a < 140°, 
the WT data is only used for a > 140°. The full CFD curve is combined directly with the WT curve for 
a > 140°. No anchoring is required since the angle-of-attack range for the WT data extends fully to 180°, 
and no modification of the CFD curve is necessary to ensure smooth blending because the two data sets are 
very close near a = 140°. Figure [IT] shows the blending for Cd and Cl • The available aerothermal data was 
not used to develop the nominal curve, but it is centered around the nominal curve. 



Figure 17: Blending IDAT CFD with 09-CA data, M = 6 
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III.F. Blending with Aerothermal Data Set, M = 8 


At M = 8, only CFD data is available. In addition to the primary Overflow set and the Usm3d set 
that cover the full angle-of- attack range, there is the aerothermal data set available in the trim range. 
Since this aerodatabase breakpoint is the border between the two database development methodologies, the 
aerothermal data set is incorporated into the nominal. This is accomplished using the same methodology 
as was used for the subsonic data, with the aerothermal data treated like the WT data and the blending 
occurring in the same manner. 

Figure p~8] shows the Overflow ID AT CFD data, the available aerothermal data, and the final nominal 
curve for Cd and Cl over the trim angle-of-attack range. The drag shows a noticeable effect of the real gas 
air chemistry used for the aerothermal data set. This effect of the real gas chemistry is incorporated into 
the uncertainty model, as described in a later section. 




Figure 18: Final nominal aerodynamics compared to available data for M — 8, a > 140°. 


IV. Uncertainty Development 

The uncertainty buildup followed the general methodology being employed by CAP. 40-45 Individual 
contributions are combined with either an RSS or additive method. The final buildup equation for the 
subsonic uncertainties is given as 


UC X — $Re T MIhs Mlidat MI (<J wt2cfd)^ T & cfd))^ T (&interp)^ i (M < 1). (2) 

The formulation changes somewhat above Mach 1 to reflect the different data used in the nominal develop- 
ment, such that 

UC X — $ltp T $ alt T $ chem T $hs T* Ml idat MI k^J '((7 wt2cfd)^ T (TTlClx((7 w i : (J cfd))^ T (&interp)^ •> (M !> 1). (3) 


The terms under the square root represent the variations seen in the data for common conditions, and 
include wind tunnel measurement /repeatability uncertainties, variations in CFD solutions, wind tunnel to 
CFD variations, and data interpolation uncertainties. These terms are combined with a root sum square 
(RSS) because they are assumed to be independent of one another (i.e. uncorrelated). They are developed 
based using range-based standard deviations that are modeled across the FMV-olt space, as described in 
detail in the companion paper 7 for the development of the cri tp term. 

Terms that are additive (e.g. Sr c ) represent a range that is covered by the database, such as the range 
of Reynolds numbers seen in flight. These terms are similar to the variation terms, but are based on ranges 
of data rather than standard deviations. 

The various margin index (MI) terms are applied to the RSS term to cover unquantified variations such 
as subsonic heatshield asymmetry M^ s , the variation due to the change in backshell angle and a 

general term covering unknown unknowns. The MI also provides coverage for lack of data. It does not cover 
large deviations from the current geometry. The development of each of the terms is covered below. 
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IV. A. Variation terms, RSS 


The terms in the uncertainty buildup that are combined with a root-sum-square, RSS , method are standard 
deviation terms for the specified variation and are expressed as, 

RSS = yj (o wt2cfd)^ T {jnax{o w t, <J cfd))^ T (o interp)^ • (4) 

These terms are developed using a range-based formula for a, 


di = Ri/d 2 , (5) 

where d 2 is a bias correction factor for converting range to standard deviation, based on sample size. 42 
The collection of cq values for a particular term are analyzed, and a model, oq, for the variation term is 
developed, typically as a function of M^, a. The model is then incorporated into the uncertainty buildup. 


IV.A.l. Wind Tunnel Accuracy , cr wt 

The term for the wind tunnel accuracy, o wt , combines a term for repeatability and reproducibility (o wtRR ) 
and balance accuracy (< o wtbal ) : 

°wt = \J X^wIrr) 2 + (vwtbai) 2 - (6) 


The balance accuracy term was developed from balance calibration data provided by the Ames Balance 
Calibration Lab 46 for the MK13A task balance used in the 05-CA test. The 1 -a error values provided in 
Ref. 46 were transformed into the force and moment accuracies using the methodology for force type balances 
described by Ulbrich. 47 The accuracies were then converted to coefficient form using the dynamic pressure 
for representative runs at each Mach number. Similar methodology was used for the 03- C A and 09- C A tests. 

The repeatability term is developed as a function of (M^a) based on the repeat run data from the 
various wind tunnel tests. For the 05-CA and 03-CA wind tunnel tests, the reported repeatability was a 
single set of values covering the entire test range. The data from these tests were re-analyzed to produce 
a model that better reflected trends, particularly in Mach number, for the test-to-test repeatability. The 
repeatability term for M = 6 was taken directly from the 09- C A report. 11 

Development of the repeatability term started with Chan’s MatLab® repeat run analysis process 48 to 
compute the individual or r. values for the cardinal angles of attack for the test at each Mach number. Only 
repeat data for the high Reynolds number data used in the nominal coefficient data was included in this 
analysis. Figure [19] shows these values plotted vs. both Mach number and angle of attack, for drag and 
pitching moment. The repeatability improves greatly for supersonic Mach numbers; this change in character 
suggests separate modeling approaches for subsonic/transonic and supersonic Mach numbers. The 05-CA 
data for above and below = 1.4 were acquired in separate test sections of the Ames UPWT facility; the 
split in the modeling approach for the repeatability is therefore made between the M = 1.4 and M = 1.6 
database breakpoints. Trends in angle-of- attack in the orr 1 data are evident for pitching moment (and lift, 
not shown) but not for drag with all of the data pooled as in Figure 

For the subsonic to low supersonic data from the 05-CA test in the 11x11 leg of AUPWT, a multivariate 
regression within MatLab® was used to fit the variations to a model that gives orr as an f(M 00: a). 
This approach allows the model to capture trends in both Mach number and angle-of-attack and provide 
repeatability data at the Mach numbers where there is not any available (e.g., = 0.6, 0.9). Figure |20| 

shows the orr model (blue diamonds) and the associated coverage, koRR , compared with the Moo = 0.50 
05-CA repeat data. The average of the ORR i data at the current Mach number is shown by the red line. 
Generally, the multivariate model provides a better description of the repeat data than the average. The 
blue line, koRR , is the coverage level for the uncertainty term, and should cover most of the data. 

For Moo > 1.4, the average of the orr 1 data at each Mach number was used. Because there were no 
significant trends with angle-of-attack, and the values for ORR aVg provided a reasonable representation of 
the ORR i data, this simpler approach was chosen. For Moo = 2.0, where there was no repeat data available, 
o RR avg for Mqo = 1.6 was used. For the Mach 6 test, the reported repeat run data 11 was used directly. 

f This term is ultimately derived using multiple levels of variation, representative of the dominate level of variation, and not 
limited to short-term measures of repeatability. 
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Figure 19: Wind tunnel repeatability range-based standard deviation data for the 05-CA and 03-CA tests. 





Figure 20: Development of model for ctrr as an f(M 00 ,a), for Co, Cl, and C m at M = 0.50. 


Sample buildups for the full wind tunnel accuracy term (Eqn. |6| are shown in Figure |2l| for Cp, Cl and 
Cm at M = 0.5 and M = 2.5. The green line with symbols gives the balance accuracy term, cr wtbar The curve 
used for cr wtRR from the previous plot set is again shown with blue diamonds in the uncertainty buildup. The 
red squares show the RSS of the two terms; this is the term that is used in the final uncertainty buildup for 
a w t. In general, the wind tunnel accuracy term is dominated by the run-to-run repeatability. 
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(a) C D , M = 0.50 


(b) C L , M = 0.50 


(C) Crr 
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(d) C D , M = 2.50 


(e) C L , M = 2.50 


(f) Cm cgx , M = 2.50 


Figure 21: Example wind tunnel accuracy buildup, for Cd- Cl and C m at two Mach numbers. 


IV. A. 2. CFD Accuracy , cr c fd 

The CFD accuracy term, Equation]?] combines a repeatability term, cr ov f RR , for the OVERFLOW computations 
with a term, cr cj g c , representing differences between Overflow and Usm3d at common conditions. 

Cc/d = yj Vox ,/ R J 2 + (o-cSc ) 2 (7) 


The set of Overflow solutions on the nominal backshell geometry had considerable duplication, with 
the differences between cases including the CFD analyst and minor code and grid inputs. While these 
solutions should be nearly identical, they were not. Their variation was used to develop a repeatability term 
for CFD in a similar manner as the wind tunnel repeatability term, 

C x -C x 

'A'max mzn 

° ov fn Ri = J 2 


There were either 2 or 3 CFD solutions available for each condition, and the repeats were only at M ^ 
= 0.3, 0.5, 0.7, and 0.9. The model for cr ov f RR was developed at these Mach numbers, and extended to 
smoothly cover the full subsonic range, up to M = 0.95. Figure [22] shows the repeat term for M = 0.5. 
The largest differences occurred in regions where the coefficients were changing more rapidly with angle of 
attack, particularly near a = 60° and a = 155°. 

The second term, cr c ^ c , was developed at M = 0.5, 1.2, 2.5, and 8.0 for the full angle-of-attack range, 
with the individual variation terms computed as 


^ OV f usmi 


Cx ov f C x 


(9) 


Some Usm3d data was available for =0.7, 0.9, and 0.95 in the trim angle-of-attack range, but was not 
used in the uncertainty development. This data had been largely discounted due to previously described 


19 ofl34l 


American Institute of Aeronautics and Astronautics 





Figure 22: Development of model for cr ov f RR as an f(a), for Cd, Cl, and C mcg at M = 0.50. 


problems with the subsonic trim CFD data in general, and because the CFD data was not used for nominal 
coefficient development in this angle-of- attack range below M = 1.6. Figure [23] shows the code-to-code 
comparison for M = 0.5. Note that the data around a = 150° is ignored in the model for a C 2 C • 


g c2c ,M=0.50 

o 

a . model 

c2c 

__k*a c2c coverage 

o 


o 



O 


°0° 

9 ~ 


% 30 60 90 120 150 180 

Angle of Attack, deg. 

(a) Cd 


0.2 


0.15 


i 

O 

b 0.1 


0.05 



Angle of Attack, deg. 

(b) C L 



(c) C mi 


Figure 23: Development of model for cr C 2 c as an f(a), for Cd, Cl, and C m at M = 0.50. 


Other options for the a c fd term were explored. Using the CFD iterative convergence criteria was not 
deemed appropriate. While the minimal convergence criteria was consistent for all CFD solutions, it was 
typically exceeded greatly for steady flow solutions. For unsteady flow solutions, where the coefficients were 
averaged over time, the convergence criteria was not a relevant parameter. Grid convergence levels were also 
difficult to use, primarily because the appropriate values were not readily available for the current CFD. 


IV. A. 3. Wind Tunnel - CFD Variations, cr wt 2 C fd 


For the subsonic Mach numbers, there is sufficient data to develop a term that compares the wind tunnel 


data to CFD data at wind tunnel conditions. This term is developed as Equation 10 


&wt2cfd i 


C Xwt Cx c/dwtRe 

d 2 


(M < 1). 


( 10 ) 


The model for this term is shown in Figure [24] for M = 0.5. Final modeling of this term (subsonically) was 
problematic due to the large jumps in the CFD data around a = 150, and the fact that this data was not 
used in the nominal development. This term was kept in the buildup, however, to provide some recognition 
that there are legitimate discrepancies between the WT and CFD data. The average value of cr wt 2 C fd . was 
used for the model, with some reduction near a = 0°, 180° to recognize symmetry conditions. This term, 
more than most, is developed with much engineering judgement. 

For the higher Mach numbers, there is not a direct comparison between the wind tunnel data and CFD 
data at wind tunnel conditions. The only available CFD data is at flight conditions, and is on the ID AT 
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Figure 24: Development of model for a wt 2cfd as an f(a), for Co, Cl, and C m at = 0.50. 


geometry. An approach is used which computes a term similar to Equation 10 using the available data, 


&wt2cfd flt . 


Cx„. - c a 


idat 


(M > 1). 


( 11 ) 


The altitude variation (which is comparable to Reynolds number variation at higher speeds) is then removed 
to provide a comparable term to the subsonic term, 


&wt2cfd — yj a wt2cfd flt a \ \lV > !)• ( 1 ^) 

This term is small compared to the magnitude of the term for M < 1. The effect of the IDAT backshell is 
still buried in the cr w t 2 cfd term. To compensate for this, is set to 1.0 for M > 1.4. For 1 < M < 1.4, 

the estimated effect of the IDAT backshell is perhaps over counted. 


IV. A. 4- Interpolation Error 

The term representing interpolation error, (Jinterp , in Equations [2j [3J and[4| was not quantified for the CM 
database but included for completeness. Development of this term is problematic. Due to the lack of data 
that would be considered representative of the flight vehicle, interpolation error would have to be determined 
for each of the constitutive data sets and propagated into the final database definition. More importantly, 
interpolation error is expected to be negligible except for a small region well outside the flight envelope. For 
a > 130°, the database point density in angle- of- attack is every 2°, and should be sufficient to capture the 
smooth trends in this region with a negligible interpolation error. Outside of this range, a < 130°, the point 
density every 5°, but the trends are still smooth over most of the domain; the exception is where there are 
steep gradients around a = 60° for the lower subsonic Mach numbers. The interpolation error is ignored in 
the uncertainty buildup, such that 


®interp 


0. 


(13) 


IV.B. Flight Condition Coverage Terms, 5 Re , 5 alu Su p , S chern 

There are several terms in the uncertainty buildup that serve to cover flight conditions that are not char- 
acterized by Mach number or velocity. These terms could have been included as additional independent 
variables in the nominal database development, but were considered small enough to be covered instead by 
the uncertainties. These coverage terms are added to the uncertainty buildup, and are developed based on 
range data. For the subsonic Mach numbers, the primary additional variable is Reynolds number; the term 
S Re is developed from CFD solutions at the bounding Reynolds number conditions. Above Mach 1, the 
data is not available to provide a Reynolds number coverage term directly. Instead, the effect of altitude 
variation,^*, developed from the aerothermal data set is used. As Mach number increases, the opportunity 
for laminar flow over the vehicle heatshield increases, and so the difference between fully laminar and fully 
turbulent flow are included as 5u p . Lastly, as Mach number increases, the effects of reacting gas chemistry on 
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the aerodynamics begin to be noticed. The reacting air chemistry is not modeled by either the Overflow 
or Usm3d flow solvers, and so coverage for this effect, 5 c h e m , is developed using the aerothermal data set. 


IV.B.l. Reynolds Number Coverage, ^R e 


The term in the uncertainty buildup for the subsonic Reynolds number increment, 5r c , was developed using 
the WT and flight Reynolds number CFD and a run-to-run variation analysis. The difference between 
computations at WT and flight Reynolds numbers was computed as 


&Rei = 



(14) 





(a) C c , M = 0.30 (b) C L , M = 0.30 (c) Cm cgx , M = 0.30 

Figure 25: Development of SR e , for Co, Cl, and C mcgx at M — 0.30. 

The computed values for SR e . at M = 0.3 are shown in Figure [25] The model for SR e is shown by the line, 
and is developed by loosely following the trends of the data. Although not shown here, the variation in 
Reynolds number decreases as Mach number increases. For drag, particularly at the lower Mach numbers, 
the Reynolds number effects are greatest near a = 0° and a = 180°. This is consistent with the flow physics, 
with the largest wake sizes at these angles and the turbulent wake being very sensitive to Rep • The lift and 
pitching moment, however, show an opposite trend with angle of attack, and the model for Sr 6 is forced to 
zero for both Cl and C mcg at a = 0°, 180° in order to preserve symmetry relations. 


IV.B.2. Altitude Coverage, S a \ t 

The altitude coverage term is developed in a similar manner as the a a i t term in the hypersonic uncertainty 
buildup in the companion paper, 7 using all of the aerothermal data for FMV < 9. The range data for each 
group of Itp aerothermal data points at a single velocity and angle of attack and computed with a single 
CFD code are used directly without dividing by the factor, such that 


Jalti 


I ^ '-E max 


(15) 


Figure 26 shows the individual S a i ti values plotted vs. Mach number. Because there are not a large number 
of data points, the model for S a i t is simply the average of all of the values. Note that the model drops to 
zero for subsonic Mach numbers. 
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IV.B.S. Laminar- Turbulent Coverage, £i tp 

All of the CFD data from the Overflow and Usm3d data sets are fully turbulent. The aerothermal 
computations are run such that there is (usually) a laminar and turbulent solution for each condition, and 
these solutions generate slightly different aerodynamic forces and moments. For the hypersonic database 
development, a variation was computed to account for the notion that the true boundary layer state isn’t 
known, and that the aerodynamics should be bounded by the laminar and fully turbulent solution^} This 
term doesn’t account for the possible effects of asymmetric transition. A similar term is computed here, 
using the range data as with the altitude term above, with 

filtPi = \Cxiam — Cxturbl- (16) 


Figure [27| shows the individual 5i tp . values plotted versus angle of attack. As with the altitude term, 
the average value is used for the final 5u p model. It is interesting to note that the highest angles of attack 
(around a = 172°) show the most variation between laminar and turbulent. As more data becomes available, 
this trend may be modeled in the future. 
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Figure 27: Development of 5i tp , for Cd, Cl, and C m( 


IV.B. 4 . Real Gas Chemistry Coverage, ^ c hem 

Above Mach 2, the vehicle can experience real gas chemistry effects in the aerodynamic forces and moments. 
The majority of the computational and wind tunnel data for M < 8 does not account for these effects. A 
limited group of computations within the aerothermal data set, using Laura, compared 5-species air and 
perfect gas air models and were used to establish an uncertainty model to account for the use of perfect gas 

g Note that the bounding assumption has not been validated. 
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chemistry in the majority of the M < 8 data. The differences for a single angle of attack at 3 velocities were 
computed as 


ficherrii ~ \^x 5sp C Xpg | . (17) 

Figure [28] shows the 5 c hemi values. In general, the real gas variation increases with Mach number, as expected. 
The uncertainty model for 5 c hem is simply the line connecting the 5 c h e mi data, anchored at zero for M = 2.0. 
For pitching moment, both the variation for C mcg and C mcgx were examined, and the maximum value used 
for the 5 c hem model. 


1-a Uncertainties for Real Gas Effects 




u(L/D) 
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Figure 28: Uncertainties to cover real gas effects where not modeled. 


IV. C. Geometry Variations 

The database nominals have been developed from a mixture of smooth, axisymmetric geometries with dif- 
ferent backshell angles. The effects of the backshell angle differences and heat shield asymmetries must be 
accounted for in the uncertainties where appropriate. The effects of other geometry variations such as the 
foot wells and windows have not been analyzed, and their effect is covered by the general margin index term, 
M. 

IV. C. 1. ID AT Differences , 

For angles of attack greater than 120°, where the database nominals were developed using data from the 
32.5° backshell geometry, the difference due the change in backshell angle to 30° must be accounted for. For 
hypersonic Mach numbers, the companion paper showed neglible differences for a > 140°; this trend should 
continue down to supersonic Mach numbers. While there were overlapping computations on both geometries 
for subsonic Mach numbers, extracting the difference in the solutions due solely to geometry proved difficult 
due to grid and CFD code execution differences. For a very limited number of cases, the difference was 
isolated, and found to be small. Because the differences were small, and only a few conditions available, a 
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multiplier approach was used to cover differences where applicable. Equation [18] gives the values for Mfl-idat 
used in the uncertainty buildup. 


Mlidat — ^ 


1.0 

if 

< 120° 

linear 

if 

120° < a < 140° 

1.1 

if 

a > 140°, M < 1.4 

1.0 

if M > 1.4 


(18) 


IV. C. 2. Heatshield Asymmetry , MI h s , ^hs 

Only limited studies, primarily for hypersonic conditions, have been done to analyze the various asymmetric 
heatshield shapes. For a select few conditions, comparisons were made with Overflow at subsonic and 
supersonic Mach numbers; the geometry differences were again co-mingled with code and grid differences, 
and so those results could not be used directly. As a result, a factor, M^, was applied to the RSS term 
in the uncertainty buildup for subsonic Mach numbers. The value of MIh s = 1.2 is based loosely on the 
Overflow results. Above Mach 1.4, the heatshield coverage term 5h s developed for the hypersonic portion 
of the database was applied. Details can be found in the companion paper. 7 


IV. D. Uncertainty Buildup for Longitudinal Aerodynamic Coefficients 

Once all of the individual uncertainty terms are developed, they are collected according to Eqns. [2] and [3j 
The coverage factor remains k = y/3 everywhere. The margin index factor, MI , is currently set as 



1.4 

if 

a < 120° 


MI = < 

linear 

if 

120° < a < 140° 

(19) 


1.2 

if 

a > 140° 



The larger margin for a < 120° reflects the fact that the database nominals are developed entirely from 
CFD data in this region. 

Figures |29f[32| show the buildup of the database uncertainty for several Mach numbers, with the legend 
for all of these plots shown in Figure 29(c)| The blue region at the bottom is the RSS term, with the light 
blue lines representing the components of the RSS term. The dark green shaded region is the Reynolds 
number coverage term that is applied for subsonic Mach numbers only. The various other coverage terms for 
M > 1 are shown in progressively lighter shades of green. The heatshield asymmetry term is in yellow, and 
represents either the 5h s term for M > 1 or the MI^s term subsonically. The ID AT coverage term is shown 
in orange, with the final MI contribution in dark red. Note that the various MI terms are formulated here 
as additive terms, so that the appropriate shaded region is ( MI X — 1) RSS. 



Figure 29: Uncertainty buildup, M = 0.5. 
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■ l DAT coverage 
Asymmetry coverage 
Altitude variation 
■Lam/turb variation 
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■Re variation 
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WT accuracy 
CFD accuracy 
A WT-CFD 

(c) Legend 


For subsonic Mach numbers, eg. M = 0.5 in Figure [29} the RSS term is dominated by the WT to CFD 
comparison term (dash-dot line), and is about half of the overall uncertainty. The Reynolds number coverage 
term (Sr Cd , green shaded) is next in importance. Note that this term increases in the trim angle-of-attack 
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region for drag. For uCl and uC mcg , the uncertainty drops as a approaches 180°, due to the symmetry 
conditions discussed earlier. Note that uCl is omitted in this figure to allow for the legend. 

All of the uncertainties drop above Mach 1, primarily because Reynolds number effects are much smaller, 
as seen for M = 1.1 in Figure |30| The WT to CFD term still dominates the RSS portion of the uncertainties. 



Angle of Attack, deg 



(a) Cd 


(b) C L 


(c) c„ 


Figure 30: Uncertainty buildup, M = 1.1, Legend in Figure 29(c) 


At M = 3, Figure [31] the effect of the heatshield asymmetry term can be seen, and is roughly the same 
size as the RSS term. The real gas chemistry effects are also beginning to contribute, primarily to the drag 
uncertainty. For the highest Mach number, shown in Figure [32j the real gas chemistry effects dominate the 
drag uncertainty. 
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Figure 31: Uncertainty buildup, M = 3.0, Legend in Figure 29(c) 
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Figure 32: Uncertainty buildup, M = 8.0, Legend in Figure 29(c) 
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Comparing the current uncertainties at M = 8 with the uncertainties from the companion paper 7 high- 
lights the inconsistencies between the two approaches. Pitching moment uncertainties are comparable, drag 
uncertainties are lower for the hypersonic approach, and the lift uncertainties are lower for the current 
approach. These inconsistencies are a topic for further development for the Orion Aerodatabase. 

IV.E. Lateral-Directional Uncertainties 

The nominal and uncertainty development presented thus far has focused on the longitudinal aerodynamic 
coefficients. In reality, there is some finite level of uncertainty in the lateral-directional coefficients (C y , 
C n , and Ci). Developing these uncertainties is a challenge, as the majority of analysis is done for (3 = 0° 
where these coefficients are either zero (for compuatational analyses), or smaller than balance uncertainties 
(in experimental data). 

The lateral-directional nominal coefficients (C y , C n , Ci) are computed directly using the usual sine/cosine 
transformations based on axisymmetric assumptions, as described in the companion paper. 7 Side-force 
uncertainty is computed as 


uCy = —uCd sin /?, (20) 

Rolling moment uncertainty was set to uCi = 0.00005 based on historical considerations, and has not been 
re-evaluated. The formulation of the yawing-moment uncertainty has evolved, and is briefly described below. 
The formulation for all of the lateral-directional uncertainties is still a work in progress, with no results shown 
here. 

The uncertainty for pitching moment was developed as a function of Mach and total angle of attack 
using Equations [2] and [3] from the previous section. The database requires that the uncertainties in pitching 
and yawing moment be provided as functions of a and (3 for the operating space of the vehicle. The 
transformations applied to uC m (aT) to get uC m (a, /3) and uC n (a, (3) have evolved over several iterations of 
the database release cycle. It has been problematic to develop a consistent formulation, primarily because 
there is not a reasonable estimate of uC n (ar ) due to data limitations. For the most recent database release, 
covering M < 8, the method for developing moment uncertainties in the full (a, (3) space was reworked to 
provide a more intuitive formulation starting with the basic pitching moment magnitude uncertainty as a 
function of aqn. Both a decrease in the basic pitching moment uncertainty and the reformulation led to a 
reduction in the yawing moment uncertainties. 

The typical sine/cosine transformation from uC m (aT) to uCm(a,[3) and uC n (a, f3) has issues at a = 
0°, 180°. The yawing moment uncertainty is identically zero here (suggesting no lateral asymmetry in the 
vehicle), and the sign of the pitching moment uncertainty changes when crossing a = 0°, 180°, both physi- 
cally unrealistic situations that led to problems with trajectory simulations. As a result, the database was 
reformulated to provide uC m (M , a, (3) and uC n (M,a,f3) directly. Developing those tables was challenging, 
and has been done inconsistently for different sections of the database. The current formulation attempts 
to provide a clear methodology to go from the (relatively) easily developed uCm(aT) to the full uC m (a, (3) 
and uC n (a,/3) by considering the relevant boundary and symmetry assumptions that can be made in the 
transformation and applying some engineering judgement to set thresholds that are not easily determined 
from the data. The uncertainties are first formulated in (aj-, 0), in order to provide an uncertainty on both 
the magnitude and direction of the total moment uncertainty. In order to formulate the final uncertainties 
as a function of (a,/?), a Monte- Carlo- like process is used to vary the magnitude and direction of the total 
moment uncertainty. This process is then applied to the current uC rn (aT) uncertainties to provide the 
desired uC m (a, (3) and uC n (a,/3 ) uncertainties for the database tables. 

V. Results 

Several aspects of the final database are explored in this section. The final nominals and uncertainties are 
compared to the experimental and computational data for selected Mach numbers to provide a sense of how 
well the database models the data. The trends in Mach number at selected angles of attack are presented to 
give a sense of how the vehicle aerodynamics change during entry. Finally, a discussion of how the database 
uncertainties contribute to trim angle-of-attack and drag dispersions is presented. These examples of the 
final database and how it is used serve to show the overall scope of the database development and reveal 
areas that can be improved upon. 
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V.A. Force and Moment Comparisons, Angle-of-attack Trends 


The final nominal aerodynamic coefficients with their associated uncertainty bands are compared to the 
available experimental and computational data for selected Mach numbers to provide a general overview of 
the final aerodatabase. Figure [33] shows the aerodynamic coefficients for M = 0.5 over the full angle-of- 
attack range. In general, the nominals were developed from the experimental data (cyan symbols) and IDAT 
computational data (orange diamonds). CFD data for the nominal backshell angle (blue symbols) for both 
wind tunnel and flight conditions as well as the Usm3d IDAT computations (purple) are shown to provide 
a reference for the code-to-code variation and Reynolds number coverage provided in the uncertainties. The 
increased uncertainty in the vicinity of a = 60° is evident for Cl and C mcg . The trim region for the same 
condition is shown in Figure 34 , for Cd and C mcg . The large variation in the CFD data is evident, with the 
uncertainties not fully covering the CFD data. This is expected based on the fact that the CFD data was 
ignored in the trim region for the nominal development, but reinforces the need for improved reliability of 
CFD methods for this type of flowfield. Note also the reduction in uncertainty for the pitching moment as 
the angle-of-attack approaches 180°. 





—Database Nominal 
■^Database Uncertainty 
o WT 05-CA data 
^ CFD, Baseline, OVF, Flight 

A CFD, Baseline, OVF, WT 

♦ CFD, IDAT, OVF 
A CFD, IDAT, USM3D 


(c) Cm cg 


(d) Legend 


Figure 33: Nominal aerodynamic coefficients with uncertainties compared to available experimental and computational 
data, for M = 0.50. 
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Figure 34: Nominal aerodynamic coefficients with uncertaii 
data, for M = 0.50, trim angle of attack region. See Figure : 



(b) Cm cg 

compared to available experimental and computational 
for the legend. 


The format of Figure [33] is repeated in Figure [35] for M = 1.2 and Figure [36] for M = 2.5. The uncertainties 
are reduced substantially for Cl and C mc3 ; Cd uncertainty remains high for M = 1.2. There is not a steep 
gradient in Cl or Cm cg at the higher Mach numbers near a = 60° as was seen for the subsonic Mach numbers. 
The uncertainty band fully encompasses the data in the trim region above M = 1.4; Figure [37] shows this 
for Cd and C mcg at M = 2.5. 



(a) Cd 




(b) ^ 
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(d) Legend 


Figure 35: Nominal aerodynamic coefficients with uncertainties compared to available experimental and computational 
data, for M = 1.2. 
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Figure 36: Nominal aerodynamic coefficients with uncertainties compared to available experimental and computational 
data, for M = 2.5. 



(a) Cd 


(b) Cm, 


Figure 37: Nominal aerodynamic coefficients with uncertainties compared to available experimental and computational 
data, for M = 1.2, a > 125°. See Figure 36(d) for legend. 


V.B. Mach Number Trends 

To show how the uncertainties change over the Mach number range, Figure [38] shows the aerodynamic 
coefficients at a = 166° plotted vs. Mach number for the full Mach number range, and again in Figure [39] 
for M < 2. It is clear that the trends in the uncertainty magnitudes are not smooth in Mach number. This 
is a consequence of developing the uncertainties independently at each Mach number rather than using some 
type of response surface, and is not a desireable feature for providing smooth trajectory simulations. This 
is an area that will be examined closer for the next release of the Orion Aerodatabase. 
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Figure 38: Nominal aerodynamic coefficients with uncertainties compared to available experimental and computational 
data, for a — 166°, M < 8 
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Figure 39: Nominal aerodynamic coefficients with uncertainties compared to available experimental and computational 
data, for a = 166°, M < 2. 


V.C. Trim Angle and Cd Variation 

At the lower Mach numbers, drag is a more important parameter than L/D for vehicle performance. In 
particular, the minimum subsonic drag is an important factor in sizing the landing parachutes. The uncer- 
tainty in drag in flight is affected by the uncertainty in Cd as well as the uncertainty in the trim angle of 
attack. In the companion paper, 7 there is a lengthy discussion on the various uncertainty parameters that 
affect trim angle of attack and lift-to-drag ratio. Figure [40 (a) shows the nominal trim angle for M < 4. The 
various symbols show the trim angle when different levels of uncertainty in pitching moment, lift, drag, and 
eg location are applied. It is important to note that, for the nominal eg location used in these calculations, 
the trim angle of attack gets as high as 175°. This is well beyond the maximum angle of attack for the 
experimental data used to develop the aerodatabase, and reflects the shift in the nominal CM eg location 
over the course of the vehicle development. The dispersed Cd is shown in Figure [40 (b)| Note that even with 
the full dispersions, the subsonic drag is roughly between 0.7 and 0.9, which is below the historical Apollo 
flight data. 
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(a) Trim angle of attack (b) Cd 

Figure 40: Trim angle of attack and Cd bounds when pitching moment, eg location, lift, and drag are dispersed. 


VI. Concluding Remarks and Future Work 

This work describes the process of developing the nominal static aerodynamic coefficients and associated 
uncertainties for the Orion Crew Module (CM) for Mach 8 and below. The database was developed from 
wind tunnel test data for the baseline geometry and computational simulations of the expanded backshell 
CM geometry, both with no asymmetries or protuberances. The database covers the full range of Reynolds 
number seen in both nominal flight and ascent abort scenarios. The basic uncertainties were developed 
as functions of Mach number and total angle of attack from variations in the primary data as well as 
computations at lower Reynolds numbers, on the baseline geometry, and using different flow solvers. The 
resulting aerodynamic database represents the best estimate of the nominal aerodynamics for the current 
CM vehicle. 

Several weakness in the current database development have been noted, and represent areas for future 
efforts that will improve characterization of the vehicle aerodynamics and reduce the uncertainties. Reynolds 
number effects for the subsonic region are not validated. This deficiency is being addressed on the exper- 
imental side with a test in Langley’s National Transonic Facility that will capture flight Reynolds number 
aerodynamics. Continuing efforts are being made to improve the computational simulations of blunt bod- 
ies with highly separated flows. The formulation of the lateral-directional uncertainties is being refined to 
provide a more rigorous physical formulation. Additional aerothermal computational solutions are being 
developed that will enable refinement of the uncertainties above Mach 2. 
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